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Abstract 

In networks the rich nodes are the subset of nodes with large numbers of links, or high 
degrees. The rich nodes and the connectivity between themselves (rich-club connectivity) tend 
to dominate the organisation of network structure. Recently there has been a considerable effort 
to characterise and model the rich-club connectivity in a variety of complex networks. In this 
paper we firstly clarify a number of terms: the rich-club coefficient quantifies the density of 
connectivity between a subset of rich nodes; the rich-club structure is the rich-club coefficient 
measured across the hierarchies of nodes; the rich-club phenomenon refers to the dynamic 
behaviour responsible for the formation of the rich-club connectivity in evolving networks; 
and the rich-club ordering discerns whether the connectivity between rich nodes in a non- 
evolving network (i.e. closed network or network snapshot) is higher than a reference network 
obtained by network randomisation. We then evaluate a recently proposed null model which is 
based on an ensemble of reference networks conserving the degree distribution of the original 
network. We remark that one should not confuse the rich-club structure of a network with 
the rich-club ordering detected by the null model. We also demonstrate that the null model 
cannot identify the dynamical mechanism that generates the rich-club connectivity. The main 
contribution of the paper is that we introduce two new ensembles of reference networks based 
on the rich-club structure for non-evolving networks. The first ensemble preserves the rich-club 
coefficient (as a function of the rank of a node) of the original network. Members of the ensemble 
exhibit similar degree distribution as the original network, and for assortative networks, similar 
assortative mixing as well. We propose that this ensemble can be used to study networks where 
assortativeness is a fundamental property, e.g. to detect the community structure in social 
networks. Analysis on the ensemble also provides a different way to interpret and model the 
evolution of social networks. The second ensemble preserves both the degree distribution and 
the rich-club coefficient (as a function of node degree). The reference networks in this ensemble 
have a similar structure as the original network. We use them to quantify the correlation profile 
between the rich nodes and pinpoint which links between the nodes are the backbone of network 
structure. 



1 Introduction 

Many social, economic, biological and technological networks contain a small set of nodes which 
have large numbers of links, the so-called rich nodes. In some networks the rich nodes are tightly 
interconnected between themselves, forming a rich-club [31] . The rich-club is an oligarchy in that 
it dominates the organisation of the whole network. In scale-free networks [1] the connectivity 
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between the rich nodes plays an important role in network functionality, for example in the trans- 
mission of rumours in social networks |i5| or the efficient delivery of data packets in the Internet |30] . 
The density of connections between the rich nodes is quantified by the rich-club coefficient [31] . 
The rich-club coefficient and its generalisations jSl HHJ EHl ESI ESI |25] have been proved to be a 
useful measure for studying complex networks. 

Many of the complex networks studied in the literature are single networks in the sense that 
each network's structure is not one of several but unique, for example the Internet. In this case we 
do not have the equivalent to a physical law to verify whether a statistical measure obtained from a 
single network is expected or extraordinary. Instead, a common technique to assess the significance 
of a property of a single network is to use the statistical randomisation method [l^ to create a 
null model. The procedure consists of using the observed network to generate an ensemble [5] of 
reference networks via randomisation. The null model is then generated from this ensemble. 

In this paper we study three ensembles of reference networks which are generated by different 
network randomisation processes. The first ensemble has been widely studied. In one recent study 
it was used to assess whether the connectivity between the rich nodes in a network is due to chance 
or an unknown mechanism. We discuss the results and clarify the confusions arising from the study. 
We then introduce the other two ensembles, which are based on the rich-club structure for non- 
evolving networks. We analyse their properties and propose useful applications for the statistical 
physics study on complex networks. Our work not only advances the understanding of the rich-club 
structure in networks but also provides new methods for studying other properties related to the 
rich-club structure. 



2 Definitions 

2.1 The rich— club coefficient and the rich-club structure 

Degree of a node, k, the number of links a node has. The rank r of a node is its position in the 
list of decreasing order of node degrees, i.e. the best-connected node is ranked as r = 1, the second 
best-connected node is r = 2 and so on. Rich nodes can be defined as nodes with large degrees or 
small ranks. The density of connections between the r richest nodes is quantified by the rich-club 
coefficient |3T] 



^r) = 4^^, (1) 

where E^r is the number of links between the r nodes and r(r — l)/2 is the maximum number of 
links that these nodes can share. If $(r) = the nodes do not share any link at all, if <I>(r) = 1 the 
nodes form a fully connected sub-graph, a clique. As a function of degree, the rich-club coefficient 
can also be given as jSj 

^ ^sfc)- <^> 

where N^^ is the number of nodes with degrees greater or equal to k and E^k is the number of 
links between the N^^ nodes. 

The rich-club coefficients <^(r) and (p{k) are related but they are not the same. The rank gives a 
unique label to each node, and the degree can be used to group nodes into subsets. If is the node 
with degree k such that -|- 1 is the rank of the node with degree k — 1, then (pi^k) = ^{r'^)- This 
is to say, {(j){k)\k = 0, . . . , /cmax} is a subset of {$(r)|r = 1, . . . , N}, where /cmax is the maximum 
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degree in the network and N is the total number of nodes. Two networks can have the same (j){k) 
and the same degree distribution P{k) for all k, but different <&(r). 

Originally the term rich-club was defined as the set of the richest nodes that are tightly inter- 
connected. This definition is vague as 'tight' is a relative concept. Recently, Valverde and Sole [27] 
proposed a criteria to define the rich-club and hence the rich nodes. The rich-club is defined by 
the existence of a crossover at kc in (f>{k) (or in ^(r)) and this crossover characterise the rich 
nodes. 

The rich-club structure refers to the density of connections across the hierarchies of nodes. It 
is given by the rich-club coefficient (f){k) for all k (or <I>(r) for all r). The rich-club structure can 
be fully defined by the degree-degree distribution, P{k, k'), the probability that an arbitrary link 
connects a node of degree k with a node of degree k' [8^ , 

^ N{k)El-lEt'^=kPik',k") 

where {k) is the average degree. Conversely, we show in the following that fixing the rich-club 
structure will constrain the degree-degree distribution |32| lllj . From the definition of the rich- 
club coefficient (j){k), the number of links that have at one end a node with degree k and at the 
other end a node with degree at least k is 

£. = E,.-E,.,. = 0W^»^-,(.+ l)^iMl(^|t±i±il. (4) 

In terms of the conditional probability P(k'\k) that a node with degree k has a neighbour with 
degree k' , 

(^max -I \ 

Y^P(k'\k)--Pik + l\k)\, (5) 
k'=k / 

where is NP{k) is the number of nodes with degree k, and the term P{k'\k) — \P{k + l\k) gives 
the proportion of these nodes that are connected to a node with degree at least k. We know that 
P{k, k') and P{k'\k) are related as P{k'\k) = {k)P{k, k') /{kP{k)) where J2k' ^i^'lk) = 1. Hence, 
fixing (j){k) for all k constrains the possible values of P{k\k') or equivalently P{k, k'). 

2.2 The rich— club phenomenon and the rich— club ordering 

The rich-club phenomenon [31] refers to the dynamic behaviour in some evolving networks where, 
if a node becomes rich, it will tend to connect with other rich nodes forming a rich-club or join an 
existing rich-club. An evolving network model can introduce such dynamic mechanisms to repro- 
duce the rich-club phenomenon and therefore generate a network with similar rich-club structure 
as the real network [71 |3l [31] [29] . One such mechanism was introduced in 1942 by Simon [26] . Si- 
mon's model was based on the addition of new nodes and the addition of new links between nodes 
that belong to the same class, where a class is the set of nodes with the same degree. Bornholdt 
and Ebel [7] have pointed out that this growth mechanism allows different growth rate for different 
classes of nodes and hence it can create a rich-club clique. Recently Krapivsky and Krioukov [TP 
showed how the inclusion of the rich-club phenomenon in evolving networks drastically constraints 
their structure. 

The term rich-club phenomenon has been mixed with the term rich-club ordering in the re- 
cent study [8] which compared the rich-club coefficient of the original network against a reference 
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network. We suggest to use the term rich-club ordering to refer to the increase of rich-club coeffi- 
cient in comparison with the reference network. It is a static property obtained for non-evolving 
networks (closed networks) or a snapshot of evolving networks. 

3 The relative rich— club reference network 

A widely studied ensemble of reference networks are the maximal random networks generated by 
randomly reshuffling link-pairs of the network under study \13\ I14j . The intrinsic structure of 
the original network is taken into account by imposing the restriction that the reshuffling process 
should not change the degree distribution. 

Recently Colizza et al. |8] used this ensemble as a null model to discern whether the connections 
between the rich nodes in a network is due to chance or due to an "organisational principle" . To 
do so, they compared the original rich-club coefflcient (l){k) with the randomised rich-club coeffl- 
cient (f)ra.n{k) obtained from the maximal random networks in the ensemble. Colizza et al. suggested 
that the "normalised" rich-club coefflcient Prwik) = ipik) / (j)ra.n{k) discounts the structural corre- 
lation imposed by the finite-size effects |6]. If /9ran(^) > 1, it means there is an organisational 
principle that leads to an increase in the density of connections between rich nodes in a more 
pronounced way than in the null model, i.e. the rich-club ordering. 

We prefer to call Pran(^) the relative rich-club coefficient. Since the number of nodes with degree 
k does not change by the randomisation procedure, pran{k) = -E^fc/-E'ran,^fc is the ratio of original- 
network links to the reference-network links. Hence Pran(fc) does not give information about the 
density of connections between rich nodes (which is measured by the rich-club coefficient). This 
simple observation is relevant because there has been some confusion in the literature on what the 
rich-club coefficient and Pran{k) are measuring. For example Fig. 1 shows the Internet network at 
the Autonomous System level (AS-Internet) O [2l] and for the scientific collaborations network 
in the area of condensed matter physics (Collaborations-A) [T71IIE]- Figs. 1(a) and (c) show the 
two networks' original (/){k) in green colour and values of the randomised (/>ran(^) obtained from 10^ 
maximal random networks, where the frequency of a particular value of (pranik) is labelled with 
different colours, from seldom (0.1, red) to often (1.0, blue). The null model (pranik)) is obtained 
by averaging over all the maximal random networks and hence corresponds to the 'bluest' dots. For 
the AS Internet, the observation that {(pranik)) > 4>{k), i.e. (pran(^)) < 1, for almost all values of k, 
suggests that the Internet does not have a rich-club ordering. This has created the misinterpretation 
as stated in [HI E] that the rich nodes in the AS Internet were not tightly interconnected with each 
other. However, Fig. 1(b) shows that the 20 best-connected nodes in the AS Internet are tightly 
interconnected between themselves. Even more, the top seven best-connected nodes form a fully 
connected clique. For the Collaborations-A network, the property {4>ran{k)) < (j){k) insinuates that 
the top scientists form tighter collaborations compared to the reference networks and this has been 
interpreted in [8l [2] as the rich nodes are tightly interconnected. However, the fact is that the top 
20 best-connected nodes are sparsely interconnected in both the original network and the reference 
networks. 

We remark that the rich-club ordering of a network is a relative property which is based on the 
comparison with reference networks generated from the network itself. It should not be confused 
with the rich-club structure which is measured and compared between different networks by the 
rich-club coefficient. Furthermore, as pointed out in [10], for both the AS-Internet and Collabo- 
rations network, the range of (pranik) increases as k increases, hence any assertion based upon the 
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relationship between {(j)ra.n{f^)) and (^{k) should be statistically tested. 

Here we show that the null model based on the maximal random networks cannot detect the rich- 
club phenomenon, i.e. whether there is a dynamic mechanism behind the formation of a rich-club. 
Consider the preferential attachment mechanism introduced by Barabasi- Albert (BA) The 
preferential attachment correlates the age of a node with its connectivity, i.e. 'rich-gets-richer'. As 
new nodes join the network the old nodes become richer. However if two old nodes do not already 
share a link, they will never acquire a new one during the network growth. This is to say, the 
BA growth mechanism is irrelevant to the formation of a rich-club. Figure 2 shows that if a BA 
network grows from a fully connected seed, i.e. a clique, it will contain a fully connected rich-club; 
if it grows from a poorly connected seed, e.g. a ring, it will have a poorly connected rich-club. The 
null model ((/>ran(^)) will produce contradictory results for the two networks generated by the same 
BA growth mechanism. The null model in this case reflects the connectivity of the initial seeds not 
the growth mechanism. 

4 The assort at ive reference network 

We can measure whether the maximal random networks used in the above null model discount 
the degree-degree correlation in the original network by examining the nearest-neighbours average 
degree of fc-degree nodes [21], knn{k) = I^^rf k'P{k'\k). If knn{k) is an increasing function of k, a 
network is assortative [19]; and if knn{k) is a decreasing function of k, a network is disassortative. 
If knn{k) is neither an increasing nor an decreasing function of k then a network is uncorrelated, 
more specifically, the degree of a node is independent of its neighbours' degree. If the maximal 
random networks are less correlated than the original network, the slope of fcnn(^) should be less 
pronounced than the original network. Figure 3 compares k^^{k) of the original network and the 
maximal random networks for four real networks: the AS Internet, the protein interaction network 
of the yeast Saccharomyces cerevisiae [13| , the Collaborations-A network, and the giant component 
of the scientific collaboration network in the area of network theory [20] (Collaborations-B). Fig. 3 
shows that the AS-Internet and the Protein networks are disassortative and their null models are 
also disassortative, whereas the two collaborations networks are assortative and their null models 
are uncorrelated networks. 

The assortative mixing is the inherent property of social networks [22 . Their maximal random 
networks do not reflect this property. In the following we define a new ensemble of reference 
networks, called assortative reference networks, which respect both the degree distribution and 
the assortative mixing of social networks. We obtain such networks by conservin £]the original 
network's rich-club structure measured by the rich-club coefficient <I>(r) as a function of rank r. 
We study the Collaborations-A network as a typical example of assortative social networks. 

4.1 Progressive rewiring 

For a given network, we start with a random network having the same number of nodes and links 
as the original network. A link is selected and rewired at randorrj^ We then evaluate the square 

^Note that conserving $(r) does not imply that the degree distribution P{k) and the rich-club coefficient (j}{k) are 
the same as in the original network. 

^We avoid creating self-loop, duplicate link, or isolated node. 
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Figure 1: The rich-club coefficient (p^k) (green) and (pranik) (red to blue) for the (a) AS-Internet 
and (c) the Collaborations network. For a given k, the range of values of (pranik), obtained from 
10^ maximal random networks, are divided into 200 bins such that both the dispersion and the 
frequency of the values can be displayed in the same graph. The frequency scale ranges from seldom 
(0.1, red) to often (1.0, blue). The null model (/9ran(^)) corresponds to the 'bluest' dots, (b) and 
(d) show the interconnections between the 20 best-connected nodes in the two networks. 
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Figure 2: The rich-club coefficient for two Barabasi-Albert networks. Both networks have 10*^ 
nodes and are grown from network seeds of 10 nodes: one seed is a ring (BA-Ring) and the other 
is a fully connected clique (BA-Clique). The colour scheme is the same as in Fig. 1. 
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Figure 3: Nearest-neighbours average degree of A;-degree nodes, k^nik), for the original network 
(green) in comparison with the average and standard deviation of k^^{k) for 10'^ maximal random 
networks (blue), (a) The AS-Internet, (b) the Protein, (c) the Collaborations-A and (d) the 
Collaborations-B networks. 
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error A = Er=i[^(^) - where $(r) is the original rich-club coefficient and $*(r) is the 

rich-club coefficient of the rewired random network. If the rewiring decreases the value of A, then 
the rewiring is accepted; otherwise it is rejected. This procedure continues until A is small. 

The assortative reference network obtained using this method not only conserve the rich-club 
structure, but also resembles the degree distribution and the assortative mixing of the original 
network (see Figures 4(a) and (b)). This is due to the structural constrain between the rich-club 
structure and the degree-degree correlation (see Eqs.|4] and |5]). Basically, preserving $(r) of the r 
richest nodes means preserving the density of connections between nodes with the r highest degrees. 
Community detection [21] is based on the comparison between the density of connections of the 
original network and reference networks. The assortative reference network could be useful for 
detecting community structures of social networks, which are inherently assortative. 



4.2 Analytical solution 

Here we give the analytical solution for the assortative reference network. From the definition of 
the rich-club coefficient <^(r) (see Eq.[T|, the number of links that the node with r shares with the 
r — 1 nodes of smaller ranks is 

^ E,. - E,.-. = ^ir) ±^ .^r-D ^'--' f-'\ (6) 

Assuming the Er links are randomly distributed between the r — 1 nodes, the probability that node 
r has a link with node r', where r' < r, is 

V{r) = (7) 
r — 1 

Thus the probability that there is a link between two nodes with ranks i and j is 

The above equation satisfies the property that X^^^ = 2L, where L is the total number 

of links of the original network. For a random network which has the same number of nodes N and 
the same rich-club structure <I>(r) as the original network, the degree of node i can be approximated 
by 

N 

h ~ "^Pij, (9) 

and the nearest-neighbours average degree of fc-degree nodes is approximated by 



i]ki=k \ J=l / 



knn{k)^— >^ \zy^P^,k, | , (10) 

where the sum with index i adds the nodes that have degree equal to k and is the number of 
fe-degree nodes. 

Figures 4(c) and (d) confirms that the analytical solution of the assortative reference network 
has similar P{k) and k^^{k) (obtained from Eqs.^and 10) as the original network. If the original 
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Figure 4: Assortative reference networks for the CoUaborations-A network, (a) (c) and (e) show 
the degree distribution P{k); and (b) (d) and (f) show the nearest-neighbours average degree of 
/c-degree nodes, k^^(k). (a) and (b) are obtained from the progressive rewiring; (c) and (d) are 
from the analytical solution; and (e) and (f) are from the rank-based preferential attachment. The 
original network is shown in green and the reference network is in blue. 
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network is assortative, i.e. high-degree nodes (with small r) tend to connect with high-degree nodes, 
we can expect from Eq. ([7| that the assortative reference network will also be assortative. However 
if the original network is disassortative, i.e. high-degree nodes tend to connect to low-degree nodes 
(with large r) , Eq. ([T]) will not capture this property as it does not favour the connectivity between 
high-degree nodes and low-degree nodes. 



Rank— based preferential attachment 



Eq. Q assumes a uniform probability for node r to attach to any node r' with r' < r. Here we 
modify Eq. ([s]) to include a preferential attachment mechanism and therefore favour the attachment 
to richer nodes with smaller r, 

^i2jV{i)/{i + l), ifi>j 
\2iV{j)/{j + l), ifi<j. 

This expression was obtained by using X^^^ i = N{N -\- l)/2 such that 2P{r){r + 1) Yl\=i ^ ~ ^r- 
This is a preferential attachment based on the relative rank of nodes. Figures 4(e) and (f) show 
that P{k) and knn{k) estimated from this rank-based preferential attachment are very close to the 
original network. 

This rank-based preferential attachment puts across a different way to study the collaborations 
network in the context of its reference network. The original network is a member of an ensemble 
which is defined by the density of connections between a referent group (rich nodes) and the total 
number of nodes. If the probability of connection between two nodes is related to their rank 
difference, then the collaborations network looks like a typical member of the ensemble. This 
provides a different way to conjecture how collaborations between scientists arise. In general, a 
scientist will prefer to work with a scientist of the same or higher status, where status is a relative 
conceplj^ 



5 The reference network preserving P{k) and 

Recently we introduced another ensemble of reference networks which are obtained by randomly 
reshuffling link-pairs of the original network with the restriction that both P{k) and 4>{k) are 
preserved for all k [32]. Consider a pair of links with end nodes ni, 77-2, and 714 with degrees 
ki < k2 < k-^ < k4, respectively. If node ni is linked to n2 and is linked to 724, we call they are 
assortatively wired. If a randomly chosen pair of links are assortatively wired they are discarded 
and a new pair is considered; otherwise the four end nodes of the links are reshuffled at random. 
The procedure is repeated for a large number of times. Reshuffling a pair of links which are not 
assortatively wired does not changes the number of links between nodes with degrees equal or larger 
than k, hence (l){k) is preserved. Notice that this randomisation procedure also conserves the degree 
distribution P{k). Fig. 5 shows that the original network and the reference networks obtained from 
this method show similar degree-degree correlation, and therefore have similar network structure. 

^This does not explain why the higher-status scientist will agree to collaborate, perhaps to carry on more work, 
or to keep his/her high status. 
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Figure 5: The reference network preserving P{k) and (j){k) for (a) the Internet, (b) the Protein 
and (c) the Collaborations-A networks. Nearest-neighbours average degree k^nik) of the original 
network is in green and that of the reference networks is in blue. For each network, 10^ reference 
networks are generated and averaged. The error bars give the standard deviation for each value. 

5.1 Correlation profile between rich nodes 

The above ensemble of reference networks allows us to assess the correlation between rich nodes in 
a network in relative to the randomised version of the network having similar network structure. 
The correlation measures how important a connection between two rich nodes is in terms of the 
organisation of the rich-club and the global structure. 

Consider aij as the zj-th entry of the adjacency matrix describing the original network, and 
a'-{m) is the ij-th entry of the adjacency matrix of the m-th reference network obtained by the 
above randomisation process, where M is the total number of the reference networks. The frequency 
probability that there is a connection between nodes i and j in the M reference networks is 

1 

771=1 

The correlation profile is obtained by evaluating bij = atj — qtj. The case hij — happens if the 
connectivity of the original network and the reference networks is the same. This case could happen 
if (1) the original network and the reference networks never have a link between nodes i and j; 
or (2) the original network and the reference networks always have a link between nodes i and j. 
The case hij — 1 means the original network has a link between nodes i and j but this link never 
appears in the reference networks. The case hij = — 1 means that in the original network there is 
not a link between nodes i and j but this link always appears in the reference networks. 

Figure 6 shows the correlation between the rich nodes in three real networks. The nodes are 
labelled by their rank. As the case hij = can represent two different situations we label the profile 
using two colour codes. The blue and white squares {hij = 0) are the links that define the backbone 
of the network. The existence, or not, of these links is fundamental for the network structure and 
perhaps network functionality 0. 
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Figure 6: Correlation profiles for the top 20 richest nodes for the Internet, Protein and 
Collaborations-A networks. For each network, the profile is obtained from 10^ reference networks 
preserving P{k) and (j){k). The colour codes are: from blue to red (6 = — > 6 = 1) labels the 
chance that a link exist in the original network but not in the reference networks; from white to 
green (6 = 0^6 = — 1) labels the chance that a link does not exist in the original network but 
exist in the reference networks. 

In the Internet profile, the 7 richest nodes are always fully interconnected with each other in 
the original network and in each of the reference networks (blue). Such a rich-club clique is an 
important structure for the Internet because it provides a large number of shortcuts for the delivery 
of data traffic and makes the average shortest path between any two nodes as small as just over 
3 hop^ We conjecture that the existence of such the rich-club clique is a fundamental property 
if the ensemble represents networks that deliver traffic efficiently. Another interesting behaviour 
in the Internet's correlation profile is that there is a link between nodes with ranks 15 and 17 
that is present in the original data but appears very rarely in its randomised version (bright red 
square). Whether this reflects erroneous or incomplete measurements, or for some reason these 
nodes interact with each other against the odds, makes this interaction of particular interest. 

The profile of the Protein network shows a completely different behaviour. The rich nodes in 
the ensemble tend not to connect with each other. The white vertical and horizontal bands show 
there are interactions between proteins that never occur in the original network or the reference 
networks. To decide if this structure is reflecting impossible protein interactions will require a more 
specific analysis. 

The profile of the Collaborations-A network shows that some researchers tend to always collab- 
orate with each other (blue) and others never collaborate (white), possibly reflecting the friendship 
and rivalries between researchers. What is never present is a collaboration between two researchers 
that it should not happen (no red squares in the profile), i.e. hindsight cannot be detected via the 
reference networks. 

*The AS-Internet contains more than 9 thousand nodes 
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6 Conclusions 



The rich nodes of a network tend to dominate the organisation of network structure so it is relevant 
to understand if their interconnectivity is due to chance or to an organisational principle. A 

technique to analyse this issue is to compare the rich-club coefficient against the equivalent structure 
evaluated from an ensemble of reference networks. The ensemble is generated via the restricted 
randomisation of the original network. 

We evaluate a recently proposed null model which is based on the ensemble of maximal random 
networks conserving the degree distribution of the original network. We remark that one should 
not confuse the rich-club structure of a network with the rich-club ordering detected by the null 
model. 

We presented a method to generate, from an assortative network, an ensemble of assortative 

reference network with the condition that the rich-club coefficient $(r) is conserved. The assortative 
reference network could be used for detecting the community structure of social networks which are 
inherently assortative. The assortative reference networks also provides a different way to explain 
the evolution of social collaborations. 

We also presented a method to generate reference networks which preserve both the degree 
distribution and the rich-club coefficient (/>(fc). The ensemble of such reference networks have a 
similar structure as the original network. We use them to study the correlations between the rich 
nodes and pinpoint which connections between them form the backbone of the network. 

Finally we'd like to remark that in the all cases of ensembles, the reference networks are gen- 
erated from the original network. Properties of the ensemble strongly depend on the properties of 
the original network. Any deduction based on these ensembles should take into consideration this 
dependance. 
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